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Response of Chestnut Flowering in Beijing to Photosynthetically 
Active Radiation Variation and Change in Recent Fifty Years 
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Abstract; Climate warming has affected plant phenology throughout the world, but few studies have evaluated plant 
phenology response to other climate factors (e. g. photosynthetically active radiation-PAR). In particular, the re- 
sponse of fruit flowering to PAR variation has not been explored yet. Long-term ( 1963-2008) of chestnut ( Castanea 
mollissima Blume) first flowering dates from Beijing, China were related with daily PAR for the 12 months, using 
Partial Least Squares (PLS) regression analysis. Two relevant phases were identified, during which mean PAR, 
temperature, and relative humidity (RH) were correlated with flowering dates, respectively. 

PAR during the both relevant periods decreased significantly in Beijing over the past 50 years. Reduced PAR 
during 24 September-5 February showed an advance impact on chestnut flowering, and could explain 12% of ad- 


vance trend in flowering timing. Deceased PAR during 6 February-31 May had a delayed effect on tree flowering, but 
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it was not significant enough to reject the null hypothesis of no impact over time. Advanced flowering of chestnut was 


mainly determined by increasing temperature between 6 February and 31 May which could explain 41% of flowering 


trend. Relative humidity variation during this period played secondly important role on tree flowering. Considering the 


interaction among these three climate factors, the impacts of PAR and RH on flowering timing could partially be at- 


tributed to the effects of temperature variation. 


Key words: Chestnut; First flowering; Temperature; Relative humidity; Photosynthetically active radiation; Partial 


Least Squares regression 


Long-term phenological observations at species 
level and specific sites can provide direct evidence 
of the impacts of climatic variation and change 
( Rosenzweig et al., 2008). While the majority of 
published studies have focused on phenological re- 
sponses to global warming and interpreted phenologi- 
cal advance as the consequence of increasing tem- 
perature (Esterlla et al., 2007; Peñuelas and Filel- 
la, 2001; Vitasse et al., 2011), relatively few re- 
ports are available on the response of plants to other 
climatic factor variations. One such change is a de- 
crease in the levels of photosynthetically active radi- 
ation (PAR) reaching the Earth’ s surface. 

PAR is the visible portion of the electromagnet- 
ic spectrum (400-700 nm) of solar radiation, cov- 
ering both energy and photon terms ( Hu et al., 
2007). PAR plays an important role in plant 
growth, development, biomass production, and 
flowering responses ( Fausey et al., 2005; Gregoriou 
et al., 2007; Hicklenton, 1987; Nathan, 1984; 
Olesen and Grevsen, 1997; Wu et al., 2010). It is 
also vital to understand the Earth’ s climate system 
and climate change (Hu et al., 2010). Advanced or 
delayed flowering events have been reported recently 
mainly through manual control irradiation amounts 
(Fausey et al., 2005; Hicklenton, 1987; Saifuddin 
et al., 2010). Responses of plant phenology to PAR 
variation in nature are reported less frequently. 

Located on the North China Plain, Beijing has 
the longest and most abundant records of phenology 
in China (Zhang et al., 2005). Compared to Ameri- 
can and European zones at similar latitude, climate 
variation and change in Beijing are substantially 
greater, leading to greater variation in the timing of 


phenological events and providing valuable informa- 


tion for elucidating climate responses of species ( Lu 
et al., 2006). Detailed and long-term meteorological 
observations in Beijing further provide an unparallel 
opportunity to assess impacts of PAR changes on 
plant phenology. Beijing currently has a total popula- 
tion near to 21 million. Rapid economic growth and 
urbanization progress have caused severe air pollu- 
tion (e. g. haze, sand and dust storm) in and a- 
round Beijing. Increased aerosol mass loading in the 
atmosphere likely affects both the quantity and quali- 
ty of solar radiation received at earth surface. De- 
creasing PAR during 1958-2005 has been reported 
in Beijing (Hu et al., 2010). How will plants re- 
spond to reduced PAR? What will happen for plant 
phenology with decreasing PAR? More scientific at- 
tention should be paid to such research to close the 
knowledge gap. 

Response of fruit tree phenology to PAR varia- 
tion, to our knowledge, has never been reported. In 
Beijing, there is a long history of chestnut cultiva- 
tion, and chestnut production generates tens of mil- 
lions of dollars in revenue for farmers. For chestnut, 
flowering phenology has vital impacts on pollination , 
fruit set and production (Legave et al., 2008). Our 
analysis was based on an observational record of 
chestnut first flowering from Beijing, China. A novel 
method ( Partial Least Squares regression) was used 
to correlate the phenology dates to PAR variation at 
daily resolution. The objective of the present study is 
to identify the relevant periods during which PAR 
can influence chestnut first flowering, and to com- 
prehensively evaluate responses of chestnut flowering 
to PAR variation during the relevant periods. The 
relative importance of other climatic factors ( temper- 


ature, relative humidity) for explaining variation in 
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chestnut first flowering has also been evaluated in 


the present study. 


1 Materials and methods 
1.1 Phenology data 

In Beijing, species-level phenological observa- 
tions of plants were mainly conducted at the Summer 
Palace (40°01’N, 116°20'E, 50m a.s.l.), a for- 
mer royal gardens with a long history of more than 
150 years. Phenological data of chestnut ( Castanea 
mollissima Blume ) at the Summer Palace during 
1963-2008 were acquired from the Chinese Pheno- 
logical Observation Network ( CPON ) -a nationwide 
system of monitoring stations that has conducted 
standardized, systematic and comprehensive pheno- 
logical observations of plants and animals across Chi- 
na since 1963. Details of the phenological observa- 
tion method have been described by Wan and Liu 
(1979) and Lu et al. (2006). In this analysis, first 
flowering phenology was used, and it is defined 
when 10% of flowers are open, corresponding to 
stage 61 on the BBCH (‘ Biologische Bundesanstalt 
Bundessortenamt und Chemische Industrie’ ) scale 
(Meier et al., 1994). 
1.2 Climate data 

Since direct photosynthetically active radiation 
(PAR) measurement has not been routinely conduc- 
ted in China, daily PAR values in Beijing during 
1963-2008 were calculated according to procedures 
used by Hu et al. (2010). Required inputs for these 
calculations were daily broadband solar radiation, 
daily extraterrestrial solar irradiance, daytime length , 
and average of the cosine of the solar zenith angles 
from sunrise to sunset. They were all from the China 
Meteorological Administration (CMA) No. 54511 
site ( Beijing Meteorological Station) which is only 
2.5 km from the Summer Palace, so that climatic 
data recorded there should closely mirror conditions 
at the observation site. To ensure the emergence of 
recognizable response patterns between PAR and 
chestnut first flowering in subsequent statistical ana- 


lyses, we subjected daily PAR to an 11-day running 


mean (Guo et al., 2013; Luedeling and Gassner, 
2012; Luedeling et al., 2013). 

To evaluate impacts of other climate factors on 
tree flowering dates, daily mean temperature and daily 
relative humidity data at Beijing Meteorological Station 
during 1963-2008 were also used in this analysis. 
1.3 Identification of relevant periods influen- 
cing chestnut first flowering 

Partial Least Squares (PLS) regression was used 
to analyze the response of chestnut first flowering to 
variation of daily PAR during all 365 days of the 
year, based on data for 1963-2008. Dependent vari- 
ables were the first flowering dates ( expressed in day 
of the year) , while independent variables were daily 
PAR values for 365 days preceding chestnut first 
flowering. PLS regression, a procedure commonly 
used in chemometrics ( Wold et al., 2001) and hy- 
perspectral remote sensing (Luedeling et al., 2009) , 
is a regression technique that can be used reliably in 
situations, where independent variables are highly 
auto-correlated and where the number of independ- 
ent variables exceeds the number of observations. 
Recent studies have proved that PLS regression can 
effectively be used to analyze relationships between 
phenology and climate variation (Guo et al., 2013; 
Luedeling and Gassner, 2012; Luedeling et al., 
2013; Ranjitkar et al., 2013; Yu et al., 2010). 

The two major outputs of PLS analysis are the 
variable importance in the projection ( VIP) and 
standardized model coefficients. The VIP values re- 
flect the importance of all independent variables for 
explaining variation in the dependent variables, with 
0.8 often used as threshold for determining impor- 
tance (Wold, 1995). The standardized model coef- 
ficients indicate the strength and direction of the im- 
pact of each variable in the PLS model ( Luedeling et 
al., 2013). Periods with VIP scores greater than 
0.8 and high absolute values of model coefficients 
represent the relevant phases influencing chestnut 
bloom timing (Guo et al., 2013). Positive coeffi- 
cients indicate that positive deviation of PAR during 


the respective period are related to delayed first flow- 
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ering, while negative coefficients imply an advanced 
impact on flowering. 

1.4 PAR temporal trend during the relevant 
periods and its impacts on chestnut flowering 

Linear regression was used to analyze the recent 
50 years’ PAR temporal trends and the relationship 
between chestnut first flowering phenology and mean 
PAR during the relevant periods confirmed above. 
Results were assessed for significance using analysis 
of variance. To investigate the relationship between 
flowering timing and mean PAR during these two rel- 
evant periods, instead of using common multiple re- 
gression, a novel three dimensional PAR response 
surface was constructed in this analysis, which was 
interpolated using the Kriging technique in the R 
package ‘fields’ (Furrer et al., 2012). 

1.5 Impacts of other climate factors during the 
relevant periods on chestnut flowering 

Temporal trends of temperature and relative hu- 
midity during the relevant periods, and the relation- 
ships between chestnut first flowering and mean tem- 
perature and relatively humidity during above phases 
were analyzed using linear regression, respectively. 
Results were evaluated for significance using analysis 
of variance. 

All analyses were implemented in the R 2. 15. 2 
programming language (R Development Core Team, 
2012). All procedures used in this study are includ- 
ed in the package ‘ chillR’ (Luedeling, 2013), 
which depends heavily on the ‘pls? package ( Mevik 
et al., 2011). 


2 Results 
2.1 Relevant periods for chestnut first flower- 
ing to variation of daily PAR 

Between 1963 and 2008, the average first flow- 
ering date of chestnut at Beijing Summer Palace was 
2 June. The 365 daily PAR values between the pre- 
vious July and June were used as independent varia- 
bles in the PLS regression, while dependent varia- 
bles were chestnut first flowering dates, expressed in 


day of the year. Based on the VIP and standardized 


model coefficients of the PLS regression, two periods 
during which PAR showed significant correlations with 
chestnut first flowering dates were identified (Fig. 1). 

The impacts of variation of daily PAR on flow- 
ering before autumn (July-September) were not con- 
sidered since a long time span from first flowering of 
chestnut (2 June) for the next year existed, and the 
response of tree flower buds to climate factors com- 
monly starts from the autumn of previous year. Be- 
tween 24 September and 5 February, model coeffi- 
cients were mostly positive and VIP values mostly 
exceeded 0.8 (the threshold for variable impor- 
tance), indicating that higher PAR should delay 
first flowering of chestnut. However, this period also 
included two spells with negative model coefficients 
indicating that the direction of PAR effect varied 
throughout this period. The lower absolute values of 
negative model coefficients meant a relatively smaller 
impact on chestnut flowering compared with the en- 
tire delayed effect of higher PAR. Taking a broader 
view at model coefficients and VIP scores during this 
period , we interpreted the entire period (24 Septem- 
ber to 5 February) as one relevant phase for the tim- 
ing of chestnut first flowering. 

Between 6 February and 31 May, most model 
coefficients were negative and VIP values were al- 
most always important, justifying consideration of the 
entire period as a homogeneous block, during which 
increasing PAR were correlated with advanced chest- 
nut flowering, while decreasing PAR linked to de- 
layed flowering. 

2.2 PAR temporal trend during the relevant pe- 
riods and its impacts on chestnut flowering 

During 1963-2008, PAR in Beijing displayed 
significant decreasing trends in the two relevant peri- 
ods confirmed above, especially in autumn and win- 
ter ( P<0. 01) (Fig. 2). Decreased mean PAR du- 
ring 24 September- 5 February advanced chestnut 
first flowering by 0.75 day per one unit decrease of 
mean PAR (P<0. 05). The correlation between mean 
PAR in later-winter and spring (6 February to 31 


May ) and chestnut flowering was not significant ( P> 
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0.1), although delayed impacts of decreased PAR 
on tree bloom were identified in Fig. 1 and Fig. 2. 
Plotting first flowering dates of chestnut as a 
function of mean PAR during the both relevant peri- 
ods clearly showed that chestnut bloom dates were 
primarily determined by variation of mean PAR in 
autumn and winter (24 September to 5 February). 
Reduced autumn and winter PAR could advance first 
flowering of chestnut. The limited influence of varia- 
tion of PAR in later-winter and spring was indicated 
by almost vertical contour lines in Fig. 2 (right part). 
2.3 Impacts of other climate factors during the 


relevant periods on chestnut flowering 





Jul Aug Sep Oct 


Dec 


In Beijing, mean temperature during the two 
relevant periods showed significant increase trends, 
by 0.5 and 0.7 °C per decade ( P<0. 01), respec- 
tively, between 1963 and 2008 (Fig.3). The ad- 
vance in blossoming was significantly related to mean 
temperatures between 6 February and 31 May. First 
flowering advanced by 2.5 days per 1 °C rise of mean 
temperature during this period ( P<0. 01). Tempera- 
tures between 24 September and 5 February were 
not significantly related to first flowering (P>0. 1). 
Rough 41% of advance trend of chestnut flowering 
can be explained by temperature variation between 6 


February and 31 May (P<0.01). 
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Fig. 1 Results of PLS regression correlating first flowering dates for chestnut at Beijing Summer Palace with 11-day running means of 


daily PAR from the previous July to June. Blue bars in the top panel indicate VIP values greater than 0. 8, the threshold for variable 


importance. In the middle panel, red color means the model coefficients are negative ( and important) , while the green color indicates 


positive (and important) relationships between flowermg and PAR. The black line in the bottom figure stands for the daily mean 


PAR, while the grey, green and red areas represent the standard deviation of PAR for each day of the year 






























































3 y y D Ay y 
528 E yr FR 5 36 48 
Period during 24 Sep-5 Feb Period during 6 Feb-31 May First flowering date 
T T 264 
> 184 © 
q S 
E 164 £ 
© 5 F 
& 14 = D 
S$ 156 > 
124 Ẹ 24 2 
~ y=-0.14x+289.19 y=-0.07x+156.53* * = E 
J R?=0.76™ "e | 184 R?=0.25" i = © 
10 on > 
T r ř r T T r r r r r r T S 
960 1970 1980 1990 2000 2010 1960 1970 1980 1990 2000 2010 Ẹ& 154 Q 
a Year Year 3 g 
F By 224 S 
È 165 165 © = 
A y=0.75x+141.82 7 y=-0.12x+155.77 5 E 
S 160] R?=0.12" =s 160] R?=001 t t 2 E= 
Om . Z E- 
© a oO ic} 
2p 1554 155 : . 2 = 
— res fa = 
Q . . a 
= < 20 = 
op 1504 150 hoe Y H iL 
| 7 PEE : g 150 
D 1454 145 
z . : = 
© . . 
= pob 140 — : f + : : ; ; : 
2 10 12 14 16 18 20 18 20 22 24 2 12 14 16 18 
ia Mean PAR between 24 Sep Mean PAR between 6 Feb and Mean PAR between 24 Sep 





and 5 Feb/mol. m2. d`! 





31 May/mol. m?. d! 





and 5 Feb/mol. m2. d`! 


Fig.2 Trends of PAR and relationships between chestnut flowering and mean PAR in the relevant periods from 1963 to 2008 at Beijing 


Summer Palace. Panels in the first and second column show trends of PAR over time and relationships between the timing of flowering and 


mean PAR during relevant periods. Trends are significant with * P<0. 1, ** P<0.05, ““ P<0. 01. The right panel displays first flowering 


dates of chestnut as a function of mean PAR during 24 September to 5 February (x-axis) and 6 February to 31 May (y-axis). Variation in 


color reflects variation in first flowering dates, while black dots indicate phenological observations between 1963 and 2008 
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Fig. 3 Trends of temperature and relative humidity (RH), and relationships between chestnut first flowering and mean temperature and 


RH in the relevant periods from 1963 to 2008 at Beijing Summer Palace. The top row is the analysis for temperature, and the bottom row for 


RH analysis. Trends are significant with “P<0.1, “ P<0.05, “* P<0.01 
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Mean relative humidity (RH) during 24 Sep- 
tember — 5 February and 6 February — 31 May has 
decreased significantly since 1963 ( P<0. 05). Neg- 
ligible impact of variation of mean RH during the for- 
mer period on tree flowering was clear in Fig. 3 (P> 
0. 1 and R’=0.01) , while decreased mean RH du- 
ring 6 February — 31 May advanced chestnut flower- 
ing significantly (P<0.01). Rough 18% of ad- 
vanced chestnut flowering can be attributed to varia- 


tion of mean RH during this period. 


3 Discussion 
3.1 Long-term trend of PAR in Beijing 

Based on an all-weather reconstruction model, 
Hu et al. (2010) evaluated the inter-annual varia- 
tions and long-term trend of PAR in Beijing. During 
1958-2005, significant decrease in PAR had been 
observed (0. 12 mol + m+ a!) which was similar to 
previous studies in Beijing (Che et al., 2005). The 
decreasing trend of PAR was sharp in spring and 
summer, but slight in autumn and winter. In present 
research, PAR during the relevant periods showed 
significant decreasing trend between 1963 and 2008, 
mainly occurring in autumn, winter and spring (24 
September — 31 May). 

In the cold winter of Beijing, large amounts of 
fossil fuels are burned for heating which has been 
exacerbated recently by rapid population and eco- 
nomic growth, and urbanization progress, causing an 
increase in aerosol emission. These aerosol particles 
can absorb and scatter solar radiation in the atmos- 
phere and likely result in the decrease of PAR in au- 
tumn and winter. In spring, the decrease of PAR 
seems to be caused by frequently occurred sand 
storms which have assailed Beijing even since Yuan 
Dynasty of China (1271-1368 A. D.) (Song, 2002). 
However, sand storms have been weakened and air 
quality has been improved since 1990 due to the im- 
provement of land surface ecological environment in 
and around Beijing (Hu et al., 2010), especially 
the construction of the Three-North Shelter Forest 
Program initiated in 1978 (Yang et al., 1998). In 


our analysis, fluctuant but whole ascending PAR in 
spring after 1990 gave certain supports to the reports 
of reduced sand storm threat to Beijing (Fig. 2). 
However, in recent ten years, haze and dust storms 
with average particulate matter diameter of 2.5 wm 
(PM2.5) have occurred frequently in Beijing and 
have aroused worldwide attention. They also reduce 
PAR in spring to some extent. To prevent threat from 
haze, sand and dust storms, China should take more 
vigorous action to reduce discharge of greenhouse 
gas, and issue a series of relative policies in afforest- 
ation, water conservation, energy restructuring , air 
pollution control, and so on. 

3.2 PAR change impacts on tree flowering 

PAR is an important climatic factor for plant 
growth, development, productivity, and flowering 
responses (Fausey et al., 2005; Gregoriou et al., 
2007; Hicklenton, 1987; Nathan, 1984; Olesen 
and Grevsen, 1997; Wu et al., 2010). In natural 
selection, plants have evolved to induce leafing and 
flowering to coincide with seasonal peaks of irradia- 
tion, maximizing PAR interception and absorption 
(De Camacaro et al., 2002; Foggo and Warrington, 
1989; Mitchell and Woodward, 1988; Williams-Lin- 
era, 2003; Wright and Van Schaik, 1994). Many 
studies have been conducted to evaluate the response 
of plant flowering to PAR variation through artificial- 
ly increased and decreased irradiation. Low PAR can 
suppress flower bud development ( Turner and E- 
wing, 1988), and cause incomplete or delayed 
plant flowering (Fausey et al., 2005; Saifuddin et 
al., 2010) , while supplemental PAR in spring can 
advance flowering of several floricultural crops ( Hick- 
lenton, 1987). Advanced flowering phase caused by 
increased PAR has also been reported in an oak for- 
est phenology study (Sierra et al., 1996). 

Limited attention is paid to explore how PAR 
variations are related to plant flowering based on di- 
rect ground observations rather than manual experi- 
ments. Long-term phenological data and detailed me- 
teorological records in Beijing provide direct evi- 


dence of the impacts of PAR variation and change on 
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tree flowering. In the present study, the relevant pe- 
riods for chestnut flowering to PAR variation were 
confirmed using PLS analysis, and it is the first time 
to report different impacts of PAR during different 
periods on tree flowering in the scientific community. 
In autumn and winter (24 September to 5 Februar- 
y), there was a positive relationship between PAR 
and tree first flowering, and reduced PAR during 
this period has significantly advanced chestnut flow- 
ering in the past 50 years ( P<0. 05). In later-winter 
and spring (6 February to 31 May) , negative corre- 
lation between PAR and flowering phenology exis- 
ted, however, decreased PAR has not showed sig- 
nificant impact on chestnut flowering ( P>0. 1). On- 
ly taking impacts of solar irradiation on tree flowering 
into consideration, chestnut bloom responses were 
relatively dominated by an advancing effect of re- 
duced PAR in autumn and winter. 
3.3 Different climate factor impacts on tree 
flowering 

Temperature and PAR are believed to be inter- 
active to determine plant flowering since increased 
irradiance could cause thermal effect to some extent 
( Berninger, 1994; Erwin, 2006; Oh et al., 2009). 
The photo-thermal ratio (PTR) model, which was 
defined as the ratio of radiant energy to thermal en- 
ergy (temperature), has been used recently to de- 
scribe the combined effects of temperature and irra- 
diance on plant growth and development ( Niu et al. , 
2000). Other climatic factors have also been tested 
to clarify which is critical to plant phenology varia- 
tion. Recently, a flowering phenology research con- 
ducted in beech and oak forests has indicated that 
the most important factors for flowering there were 
soil moisture, and relative humidity, followed by 
temperature, numbers of daylight hours, and PAR 
(Sierra et al., 1996). However, most of current 
knowledge about plant phenology, including numer- 
ous experimental studies ( Menzel and Fabian, 
1999; Price and Waser, 1998) , indicates that ob- 
served phenology changes are mostly due to increas- 


ing temperatures (Borchert et al., 2005; Esterlla et 


al., 2007; Pefiuelas and Filella, 2001; Zheng et 
al., 2006). 

In our analysis, rough 41% of advance trend of 
chestnut flowering can be explained by temperature 
variation during 6 February - 31 May (P<0.01), 
while decreased relative humidity during this period 
can explain 18% of tree flowering trend. Comparative- 
ly, only 12% of flowering trend can be attributed to 
variation of PAR during 24 September — 5 February. 
Thus, our results were in line with most reports that 
plant phenology are mostly determined by temperature. 

There seems to be some correlation and interac- 
tion among temperature, PAR, and relative humidi- 
ty, while the latter two can both be explained 
through temperature effects. Based on long-term re- 
cords of chestnut first flowering and daily tempera- 
ture data between 1963 and 2008 in Beijing, Guo et 
al. (2013) confirmed the relevant periods for flo- 
wering to daily temperature variation which were al- 
most identical to these periods in this analysis. This 
indicated similar impacts of temperature and PAR on 
tree flowering. Decreased PAR in autumn and winter 
seems to function as reduced temperature, while the 
latter benefits the chilling accumulation of plants, 
and leads to earlier dormancy breaking and flower 
sprouting. The advanced impact on flowering of de- 
creasing PAR in autumn and winter actually oc- 
curred in this analysis, probably due to general cor- 
relation between radiation and warmth. Between 6 
February and 31 May, decreased relative humidity 
can link to increased temperature in the same period. 
Thus, the advance effect of relative humidity varia- 
tion in spring could also be explained by the corre- 


sponding temperature increase. 


4 Conclusions 

Response of fruit flowering to PAR variation has 
been evaluated in this study for the first time, based 
on long-term phenological and meteorological records 
in Beijing, China. Two relevant periods were identi- 
fied during which PAR had completely opposite im- 


pacts on tree flowering. In autumn and winter (24 
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September to 5 February), observed decreasing 
PAR advanced chestnut flowering significantly, 
while reduced PAR in later-winter and spring (6 
February to 31 May) imposed negligible effect on 
flowering. Comparatively, the most important climate 
factor for chestnut flowering was temperature, fol- 
lowed by relatively humidity and PAR. Impacts of 
PAR and relative humidity on tree flowering could 


partially be explained by temperature effects. 
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